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Abstract 



O Simulation results for Mobile Ad-Hoc Networks (MANETs) are fundamentally governed by the underlying 

Mobility Model. Thus it is imperative to find whether events functionally dependent on the mobility model 

^ 'converge' to well defined functions or constants. This shall ensure the long-run consistency among simula- 

tion performed by disparate parties. This paper reviews a work on the discrete Random Waypoint Mobility 
Model (RWMM), addressing its long run stochastic stability. It is proved that each model in the targeted 

I I discrete class of the RWMM satisfies Birkhoff 's pointwise ergodic theorem |[T3l . and hence time averaged 

l-H 

^ functions on the mobility model surely converge. We also simulate the most common and general version 

t/3 of the RWMM to give insight into its working. 
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Q^ Introduction 

en 
en 

Q Mobility models are used for the generation of node movement in simulations of MANETs. Protocol 

^ development is a consequence of such a simulation. The probabilistic aspects of the founding mobility 

^ model has direct implications on the simulation results. Many papers Q-im have already concluded that 

k> stochastically unstable mobility models shall result in simulation results that diverge in time. 

The Random Trip Mobility Model, through the presence of a unique stationary distribution for the 
location of nodes, has already been proved to be stable [5|. The work presented in this paper is purely a 
review of the stability of the discrete version of the RWMM proved by Timo, Blackmore and Hanlen |[T1. 
Therein the notion of stability is considered to be the satisfaction of Birkhoff 's Pointwise Ergodic Theorem 
by the mobility model. If to the contrary the mobility model is unstable, the simulation results are bound to 
be unreliable. 

The stimulus for this line of work is that the stability or lack thereof of the mobility model is possibly 
passed up the layers of the protocol stack. For instance the DSR protocol preserves the mobility model's 
stability IS : if the node location random process is stable, then so is the route selection random process. 
The consequence of this is that the strong law of large numbers also holds for the simulations at the network 



layer. 

A mobility model is quantified using a random process. It is stationary if the set of probability laws 
regulating the movement of the nodes are independent of time. Many works have come up with the trans- 
formation of non-stationary models into (in some places pointwise ergodic theorem satisfying) stationary 
models with the motivation that the strong law of large numbers may be applicable. 

The classic RWM model does display starting transients and local nonstationarity. Thus we analyze its 
properties by means of imposing a mathematically weaker 'asymptotic stationarity' property. A random 
process, the mean of which is governed asymptotically by a process with a stationary distribution is called 
Asymptotically Mean Stationary (AMS). It has been proved HI [Theorem 1] that a random process is AMS 
if and only if it satisfies Birkhoff 's pointwise ergodic theorem. By consequence a mobility model is stable 
if and only if it is AMS Q. 

In the classic RWMM m, every node, using an independent and identically uniformly distributed 
(IID) random process {VF^J^q, selects a sequence of waypoints w = wo,wi,W2, ■■■■ For every pair 
{'Wi,Wi^i),i G Z*, the node chooses a speed randomly and uniformly distributed from the closed interval 
[7nin_s, max.s]. At this chosen speed it then travels in a straight line from Wi to Wi+i. 

In this review, the main result addressed is: a) The general discrete class of the RWMM is asymptotically 
mean stationary (by virtue of which it is stable) and ergodic. b) For stable node movement the following 
conditions suffice - (i) Node waypoint selection is an AMS random process, (ii) Speed selection random 
process is stationary. 

This paper is organized as follows. The next section introduces the preliminaries. Following that we 
describe the general RWMM. Next is the contribution of the base paper in the form of a theorem. Simulation 
results for the classic case and conclusion end this paper. 



Preliminaries 

We will adopt the dynamical systems ||9]|- ||T0| model for a random process. Given a discrete finite alphabet 
X, let X = {Xfej^Q be the associated discrete time random process. The distribution of {X„}^q is the 
set {/i(^) : k > 0} where /x^*^) is the probability measure on X'^ given by: 

^W(xg-i) = Pr[Xo = xo,Xi = xi,...,Xk^i = Xk-i] 

In order to simplify our work, we use the Kolmogorov Representation Theorem (where certain consis- 
tency conditions are satisfied) |9] [Theorem 1.1.2]. This enables us to replace the distribution with a unique 
probability measure fi on the space X°° = U^qX. Throughout we shall be dealing with cylinder sets as 
elementary events: [x^] = {5L : Xi = Xi,m < i < n}. The cr-algebra Tx°° is generated using these 
cylinder sets. Time is incorporated using the shift transform T^ = Xk, Xk+i, Xk+2, ■■■, k G Z*. Eventually 



we result with the dynamical system {X'^,Tx°° , M, T^) which is related to the original random process by 

{Xk}r=o = {no(r|x)}^^o, nox = xq. 

Suppose we have a mobility model quantified by the random process {X^j^q and capture the location 
of each node for the first k time instances of a simulation given by Xq^^ = xq, xi, X2, ..., x^-i. The dy- 
namical system associated with this stochastic experiment is {X'^jTx"^ , M, T^) and the trajectory captured 
is the elementary event [xq~^] G J^x°^- If variable length shift must be considered Tx*x = T^ x as 
is necessitated in certain cases by the random processes associated with the updation of routing tables of 
network routers, we may study the probabilistic properties of {X°° , Tx°° , fJ-, Tx*)- 

Now we come up with certain definitions and lemmas lifted from the base work which serve as the 
foundation for future proof developments. 

Definition 1 (Stationarity)^: The system {X°° , Tx°° , /x, T^t) is called stationary, and Tx is said to be mea- 
sure preserving if, V^ G Tx°°, m(^) = /^(^ "'^^) ■ 

Definition 2 (Ereodicitv } ifTl lfT3l : The stationary system {X°°,Tx°°, ^J-jTx) is said to be ergodic if ^ = 
T-^A^ H{A) = 0or 1 . Equivalently, (A""", J'^'-', /i,7>) is ergodic iff V/ G ^^(^u) the limit 

n-l 

< / >=< / > (x) = lim - J2 fi^x^) is a constant almost everywhere in ^. 

n-s-oo " ^^Q 

Definition 3 ( Stability )f[\: A mobility model associated with the random process {X°°,J^x°°,fJ-, Tx) is said 
to be stable if for all bounded and measurable /, the limit 

n-l 

< / > (x) = lim - J2 f{Tx'^) exists almost everywhere in /z. 

Definition 4 (Asymptotic Mean Stationarity)[l\: The system {X°°,Tx^, fJ-,Tx) is said to be asymptotic 

n-l 

mean stationary (AMS) if, V^ G Tx°^ the hmit Jl{A) = lim ^ X) K^x^) exists. 

Here the probability measure Jl is defined on the measurable space {X°°,Tx°°)- It is called the stationary 

mean of fi and describes the average of the long run behaviour of the system. 

Lemma 1 (Birkhojf's Pointwise Ergodic Theorem) |[T| |[T3]| : Let the dynamical system {X°°,Tx^ ,fi,Tx) 

have Tx as a measure preserving map, and let / be measurable with E{\f\) < +oo. Then 

n-l 

lim - J2 f{Tx^) = E{f\C). Here C is the cr-algebra of invariant sets of Tx ■ If the random process is 

n-s>oo " ^^Q 

ergodic, then C is the trivial c-algebra, and E{f\C) = E{f) which is a constant. 

It has been proved JSl that asymptotic mean stationarity is both a necessary and sufficient condition for the 

pointwise ergodic theorem. 

Lemma 2 (AMS Pointwise Ergodic Theorem) jSl [Theorem 1]: A dynamical system {X°° ,Fx°^ , Z^, Tx) is 

AMS iff for all measurable / with a finite expectation, the limit 

n-l 

< / > (x) = lim - Yl, f{Tx'^) exists almost everywhere in fi. 
Eventually we conclude, using definitions 3,4 and lemma 2: 

Stability: A mobility model with {X°° , Tx-^ , /u) as the associated probability space is stable with respect 

to Tx iff (A"^, J"a- , M, Tx) is AMS. 



Discrete Version of the RWMM 

We now initiate the study of a discrete time space version of the Random Waypoint MobiUty Model. Con- 
sider a MANET with each mobile node in the set V = {vi ,V2, ■■■, ^ivi } located in a discrete finite geograph- 
ical area described by the set S. The following are the random processes to exposit the discrete RWMM. 

Waypoint Random Process Per Node 

From the geographical space S, each mobile node v G V selects an infinite tuple of waypoints w = 
wo,wi,W2, ... randomly. Let us denote the waypoint selection random process as W = {VFfcj^Q with the 
corresponding dynamical system as (>V°°, J^vv^o , ^^, T\v), and also W = S. 

RWMM Correlation : In the classic RWM model, waypoint selection random process is IID, and in most 
cases uniformly distributed. So the stochastic process (W°°, -Fyyoo, /iw, Tyy) is a Bernoulli Scheme |[T4l . 

Path Random Process Per Node 




Pi 
p. 



Figure 1 : Different paths corresponding to discrete time-space equivalent of different speeds 



In the classic RWMM, whenever an arbitrary node selects a sequence of waypoints w, then for each 



consecutive pair {wi,Wi+i),i S Z*, it also selects a speed uniformly distributed from [min_s,max_s] 
and traverses the straight line path between Wi and tUj+i. In discretized time and space, snapshot of the 
node's position per instance of time is taken during its trip between waypoints. This shall result in a random 
path with finite possibilities per waypoint pair {wi,Wi+i) (figure 1). For each combination of waypoints 
{w, w') ^ W y^yV construct the set of all paths Vw^w' and take the union of all such sets so as to obtain 
all admissible paths P: 7^«,,«,' = {pi-,P2-,P3-, ■■■,P\v ,\}, V = U V^.w- In order to de- 

scribe the stochastic process P = {-Pfcj^o' noting that P is conditioned on W, we first define the set of 

permitted path sequences V^ C P°° given vi^ as P^ = {p G V^ : Pi G 'Pw„w,+i,^'i G 2*}. Here 

again let Tpao be the a-algebra generated from p^ G V°°. Defining a collection of probability measures 

zy^p = {zy^ : w G W°°, i/wl'Pw ) = 1} ^esults in the channel iHJ (W, iy^p,'P). 

Definition 5 (Stationary Channel){\\: If Vw G >V°°,Vyl G F-poo, vty^^[A) = v^{T^'^A) , the channel 

(W, i^wp, 'P) is said to be {Tw;, T-p) stationary. 

RWMM Correlation: Considering elementary events [Po ^] G F-p^: 

n—l -■ 

^w([pr']) = n w 1 Pi e ^-.,-,+1,0 < i < n - ]l 

And Vyf is zero otherwise. To prove stationarity, see that on the transformed r>vw the non-zero probability 
for \pq'^] is given by: 

n— 1 1 
2^Tww([Po"^]) = n Vb 1 Pi ^ Vn^.+uw.+i^O <i<n-l 

And if Tp ^ [p^~^] = \p^] with p^+i =Pm,0<m<n: 

n—l 1 

^w(Tp^b^^]) = u^m]) = n p^ 1 P^ e ^«,,+„«,,+2, <i<n-i 

j=0 I Wi+i,Wi+2\ 

The last two equations are equal which proves that the channel is stationary. 

Next it is proved that the channel is output mixing and consequently ergodic. A channel is said to be 
output mixingQ if, VA, B G FpoJQ Vw G W°°: 



lim li^w fe^^fl^) - i/w(Tp"A)i/w(S 



The elementary events in case of the general RWMM are decoupled for r > 6 for [p^ ^], [p^ ^] G Fpoc in 
the following equation: 




6-1- 



p;+"-'] nipo 



Pi G 'Pw,+-i_,w,+2 ,T <i<T + a-l 
Pj e P^j.+i,^j+2, <j<b-l 



'Product upper limit error in fT| 
^Incorrect sigma field in 1 1 1 



Hence the channel is output mixing and ergodic ITTI [Lemma 9.4.3]. 

Finally we define a probability measure ^p conditioning it on the waypoint selection probability mea- 
sure fi^: fip{A) = J2 Aiw(w')i/w'(^)) ^^ £ J-poo . Thus we result with {V°° , J^-poc , fipjT-p) as the 
corresponding dynamical system for P = {Pfcj^Q. 



Location Random Process per node 

We define the time taken t{i) to reach Wi from wq as a function of the first i paths po,pi, ...,pj_i. We 
assume that each path length l{p) is a positive finite quantity. So t{i) = J2l^o KPj)^ ^ > 1- Let the i^^ 
path Pi take the form St(i), Si(i)+i, 54(^)4.2, ...,St(j+i), with Sj G S and Sj(fc) = w^. Correlating with the 
given paths' sequence p, we arrive at node location sequence s = so,si, .... Thus we have the node location 
random process S = {Sk}'kLo given by the dynamical system (5°° , J^s°° ,l^s,Ts). 

Location Random Process for all nodes 

Consider the |V| tuple X„ = (S'„,^i, 5„,2, •••7 "^n ,|v|)' '^ith Sij denoting node j's location at time i. This 
random variable's alphabet is A" = S'^K Define {X'^,J^x°°, fJ'p,Tx) as the dynamical system for the 
random process {X^j^q. 

Main Result and its Proof 

Theorem IH: Suppose the nodes V = {^1,^2, •••)^|v|} move in agreement with the discrete RWMM 
already defined. Let W^, denote the waypoint selection random process for node v and P^, be the cor- 
responding path random process. Let (Wv, i^wp,'Pv) denote the path and waypoint stochastic processes' 
connecting channel and let X denote the location random process for all nodes. Then 

• Uyv £ V, Wy is AMS and the channel (Wv, J^wp, 'Pv) is stationary, then X is AMS and stable. 

• If Vw G V, Wt, is ergodic and the channel (Wv, z^wp, Vv) is ergodic, then X is ergodic. 

Proof Sketch Q: 

Dropping the redundant subscript v henceforth. 

Lemma A: IfW is AMS and ergodic and (W, i^-wpjV) is stationary and ergodic, then P is AMS and 
ergodic. 



Proof: mi [Lemmas 9.3.1, 9.3.3] prove this lemma directly as the AMS and ergodic waypoint random pro- 
cess and the path random process are connected by a stationary, ergodic channel. 

Lemma B: IfP is ergodic then S is ergodic. 
Proof: Given (P°°, J'poo , /Xp, T-p) is AMS. For aWp £ V let l{p) denote path length, let L = max{l{p) : 
p € V} and let / : "P — )• Uj=i '5* be the breakdown of a path to its corresponding geographic cells - f{p) = 
so,si,S2,...,Si(p)_i. So S = {5fc}^o = f{Po),f{Pi),f{P2),--- = 5'o,5'i,...,5/(P(,),...,S'i(P(,)+i(Pj),.... 
For ease of working, define the encoder F : V°° — )■ 5°° as s = -F(p) = f{po), /(Pi)) f{P2), •••■ Now 
\/A G Ts°° ) ^s(^) = fJ-p{F'^A)- Here the mapping F is many to one. 

Next it is described a pseudo-inverse G~^ : S°° — )• V°° as G^^s = ps where V°° C 7^°° according to 
im [Theorem 1] having full measure //p('P°°) = 1 such that every bounded measurable function on this 
set converges; 5°° is the induced range of F on V'^ and ps is a representative from the partition of V'^ 
induced by s G S°° . 

Define the length of the first n paths in p G V°° as 7p(n) = J2^=o KPiVl Then the variable length shift 
Ts* : 5°° ^ 5°° is given by T^,s = tJ"^^^s where 

r.(s)= ^o-Hs)i^), ses* 

1, s^S* 

Eventually it is proved that lim„_i.oo ^ J^kZo h{Ts*s) exists Vs G S* and for all bounded measurable h. 

Thus (5°° , J^s°° > Ms , Ts* ) is AMS. Note that one Ts* shift is equivalent to one path shift. 

Sublemma: If (5°° , J^s°° ; /"si Tg* ) w AM^ w/?/j stationary mean /x* ?/ie?i (5°° , J-5°° i A's, TI5) /s AM^ where 

Ts'S = Tj^^^s otJ 1 < 7(s) < L . 15 

Define a new measure (inspired from HI [Ex. 6]): 

Here A^^ = {s G 5°° : 7(3) = k}, and {A^^}^^^ is a partition of 5°° [S]. Thus we have Tg}A = 
UJ^^^iT^^A n A^i). We also have 

KiTs'A) = fiHA) = 5: /z:(A n A^i) = 5: a:(T5 '^^ n A, 1) 

fc=i fc=i 

The first two terms are equal by virtue of transformation invariance of /i* . The next two terms are equal by 
virtue of intersection distribution of A on A^^. The first and the last term are equal from the immediately 
preceding correlation between Tg* and Tg. Substituting T^^A for A in the definition of /Xg and using the 
above equation we arrive at the T5 invariance of /Xg- Further it is shown that fig asymptotically dominates 
j2l under T5 which when taken with the T5 invariance of p,s and [8j [Theorem 2] proves that /i* is AMS 



'Wrong limit in fTI 
^Typographical error in 1 1] for Ts' 



w.r.t. Ts . 

Next it is proved that if /is(^) = and T^^A = A then /is(^) = 0. This in conjunction with |[T2l TTheorem 

2.2] proves that ^Ug is AMS w.r.t. Tg. 

Thus from the sublemma (5°°, J^s°° > /^s, T!s) is AMS which completes the proof. 

Lemma C: IfF is ergodic, then S is ergodic 
Proof: Let A G 7\s°° be Tg invariant and let s = F{p) be an arbitrary member of A. Now 
F(p) £A^ TgF{p) £A^ F(p) £A^ T^^^" V(p) €A^ F{p) €A^ F{Tpp) G A 
Taking F^^ on both sides (as the equation holds for all s and all p associated with each s ) 
p E F~^A <^ T-pp e F~^A . Hence F^^A is T-p invariant. By the premise of the lemma, fip{F~'^A) = 
or 1. Hence /is(^) = or 1. Thus by definition, {S°°,Fscxi , /ig, Ts) is ergodic. 

Simulation 



We have simulated a basic packet exchange in a MANET using NS2 with the node movement generated 
according to the general continuous RWMM. The traffic consisted of constant bitrate UDP packets with 
IEEE 802.11 protocol at the MAC layer. The exchanges resulted in bursty traffic. The plot for the number 
of bytes received as a function of time for a particular node is given in figure 2. 
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Figure 2: Number of bytes received as a function of time for a particular node 



Conclusion 

In this paper we have successfully demonstrated that the discrete general RWMM is AMS and hence stable. 
Thus simulations with RWMM as the underlying node movement generation algorithm tend to be reliable. 
The stability preserving protocols allow higher layers of the protocol stack to propagate this stability hence 
permitting reliability of simulations at higher levels also. Moreover we have simulated the continuous 
version of the RWMM with the intent of seeing the local non-stationary properties (which is highlighted by 
the bursty traffic). 
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